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Abstract 

Blue primordial power spectra have spectral index n > 1 and arise naturally 
in the recently proposed hybrid inflationary scenario. An observational upper 
limit on n is derived by normalizing the spectrum at the quadrupole scale and 
considering the possible overproduction of Planck mass relics formed in the final 
stage of primordial black hole evaporation. In the inflationary Universe with 
the maximum reheating temperature compatible with the observed quadrupole 
anisotropy, the upper limit is n = 1.4, but it is slightly weaker for lower re- 
heat temperatures. This limit applies over 57 decades of mass and is therefore 
insensitive to cosmic variance and any gravitational wave contribution to the 
quadrupole anisotropy. It is also independent of the dark matter content of the 
Universe and therefore the bias parameter. In some circumstances, there may 
be an extended dust-like phase between the end of inflation and reheating. In 
this case, primordial black holes form more abundantly and the upper limit is 
n = 1.3. 
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1 Introduction 



Is is now generally accepted that large-scale structure in the Universe arose from the 
growth of small density fluctuations through gravitational instability. A determina- 
tion of the power spectrum of these primordial fluctuations is therefore of fundamental 
importance and the inflationary scenario provides an attractive, causal mechanism 
for producing such a spectrum (see for a recent review). During inflation the Uni- 
verse is dominated by the self-interaction potential V{(f)) of a quantum scalar field 
and this results in a superluminal expansion of the scale factor. Quantum fluctua- 
tions in this field are therefore stretched beyond the Hubble radius H~^, where they 
remain frozen until reentry during the radiation- or matter-dominated eras. 

In general, when scalar density fluctuations reenter the Hubble radius, their am- 
phtude is given by 

1 1/3/2 

(1-1) 

where the quantities on the right-hand-side are evaluated when the scale first crossed 
the Hubble radius during inflation, a prime denotes differentiation with respect to 
and Mpi is the Planck mass [Q. The form of these fluctuations can be parameterized 
by the spectral index n which specifies the dependence of the power spectrum on the 
comoving wave-number k (viz. |(5fc|^ cx fc"). This means that the density perturbations 
have an rms amplitude of the form 5{M) oc M*^^~"^/^ when the mass scale M reenters 
the horizon after inflation. 

In the simplest scenarios these horizon-scale fluctuations are expected to be almost 
scale- invariant with n ^ 1 However, observations of large scale structure suggest 
that scale-invariant fluctuations may not work, at least for the cold dark matter 
(CDM) model 0, ||, and it may be necessary to consider 'tilted' spectra with n ^ 1 
'Blue' primordial power spectra have n > 1 and are currently consistent with 
the recent anisotropy measurements of the Cosmic Microwave Background (CMB) 
radiation [6-13]. A best-fit of the theoretical and observed autocorrelation functions 
with the COBE/DMR first-year maps implies n = l.lbtoH 0- Smoot et al. 
and Torres [§ have independently performed a topological analysis of this data and 
deduce n = l.Tlg^ and n = 1.2loi, respectively. On the other hand. Bond 
infers n = l.Stoi from ^IRS data at 3.8° ^ and n = 2.0l:oi from the first-year 
COBE/DMR data. The second-year COBE/DMR. data has now been analyzed by 
Wright et al. ^ who find n = lAGtHl, while a maximum-likelihood analysis of 
this data by Bennett et al. |jl2| implies that n = 1.59lo;55. The 53 and 90 GHz 
second-year data has also been analyzed by Gorski et al. [0, who find a maximum 
likelihood value of n = 1.22 (1.02) if the quadrupole is included (excluded). They 
further obtain a marginal probability distribution with a mean of n = 1.10 ± 0.32 
(0.87 ± O.36).0 It is interesting to note that n > 1 is deduced by Piran et al. from 
large-scale structure considerations |]T^. If the voids detected in the CfA survey |T5| 

^It is now recognized, however, that the Doppler peak influences the low multipole harmonics of 
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on scales 50h ^ Mpc arise from an underdensity in the matter distribution and form 
gravitationally in an = 1 Universe, a spectral index of n ~ 1.25 is consistent with 
the COBE/DMR detection. 

The above results implicitly assume that the spectral index is constant over the 
scales of interest, but in general n is a function of scale and is determined by the 
magnitude of the potential and its first and second derivatives fl^. However, scales 



relevant for large-scale structure lie in the range 1 - 10^ Mpc and correspond to only 9 
e-foldings of inflationary expansion. Since the scalar field must roll down its potential 
very slowly during inflation, only a very narrow region of V is relevant for these 
scales. It is therefore consistent to expand the potential as a Taylor series about a 
given scale and this is equivalent to assuming that the spectral index is constant over 
a sufficiently small range |18[. 

The full class of general potentials leading to spectra with constant spectral index 
has now been derived |19|, The potential leading to n > 1 is a combination of 



trigonometric functions [|l^ and its Taylor expansion to quadratic order is given by 



k2 



l + 27r(r2-l) 



Mil 



(1.2) 



where Vq is a constant that can be normalized to the COBE/DMR quadrupole via Eq. 
( p.. 11 ). It should be emphasized that any potential exhibiting a Taylor expansion of 
this form over scales corresponding to large-scale structure leads to a blue spectrum. 
Blue spectra also arise in some variants of the hyperextended scenario ||2T|| . 



Such potentials form the basis of the hybrid inflationary scenario |22] and can 



arise in string physics [^. The COBE/DMR detection implies that the last stages 



of inflation must have occurred at or below the Grand Unification scale and the 
superstring is an effective = 1 supergravity theory at these scales. It can be shown 
that under fairly generic circumstances the one-loop Kahler potential derived from 
the orbifold compactification of the superstring is given approximately by Eq. (|1.2| ) 
where the spectral index is determined by fundamental string parameters: 



The dilaton field S is the real component of the chiral multiplet and its vacuum ex- 
pectation value determines the string coupling constant (Re(S')) ~ S'gut • Results 
from LEP imply that (^gut ~ 27r/13 and consequently (Re(5')) ~ 2 ||2^. The pa- 
rameter S^^ is a dimensionless coefficient to the one-loop corrections arising in the 
Green-Schwarz mechanism and calculations suggest that df^ < 47r^ These values 
lead to the upper limit n < 1.5, although the requirement that inflation ends at the 
correct time leads to a stronger limit n < 1.3 

the temperature anisotropy in such a way that the observed spectral index is tilted slightly to the 
blue end, i.e. the spectral index as measured by COBE/DMR is not the index of the primordial 
spectrum. For standard CDM the correction is of order 15%. 
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We therefore infer that values of n > 1 are not excluded at this stage, either 
observationally or theoretically, and it is therefore important to obtain upper limits 
on the spectral index. Blue spectra introduce more short-scale power and this might 
be problematic for dark matter models of galaxy formation. However, non-linear 
effects become more important on larger scales when n > 1 and these effects must be 
adequately accounted for before such spectra can be excluded by large-scale structure 
arguments. In principle, the extra short-scale power can be significantly reduced by 
the free streaming of a hot dark matter component and comparison of the mixed dark 
matter model with observations of large-scale structure above 1 Mpc implies an upper 
limit of n < 1.35 ||2^, with dependence on the current value of the Hubble constant. 
The limit on the spectral distortion of the CMB from the COBE FIRAS experiment 
leads to a slightly weaker limit of n < 1.56 |^ . 

One of the most interesting constraints on inflationary scenarios which produce 
blue spectra comes from considering the formation and subsequent evaporation of 
primordial black holes (PBHs) |T^. PBHs are never produced in sufficient numbers to 
be interesting if n < 1, but they could be if n > 1. Such limits are interesting because 
they constrain the power spectrum, and therefore the inflationary potential, over 47 
decades of mass, whereas large-scale structure measurements span only 10 decades. 
This implies that the limits on the spectral index derived from PBHs are extremely 
insensitive to the precise value of the quadrupole anisotropy and so the problems 
associated with cosmic variance on these ultra-large scales are evaded. They are also 
independent of any gravitational wave contribution to the CMB anisotropy and, when 
the spectrum is normalized to the COBE detection, they are also independent of the 
transfer function, the bias parameter and the value of Hubble's constant. Although 
such observational limits are essentially independent of observations, they are derived 
from the assumption that the spectral index is constant over a very large range of 
scales. In the Appendi x it is shown that this is a valid assumption in the hybrid 
inflationary scenario. 

In this work we shall consider the potentially stronger observational constraints 
which arise if evaporating PBHs leave stable Planck mass relics. Several people have 
considered the cosmological consequences of such relics. MacGibbon considered the 
possibility that they could have around the critical density and thus provide the dark 
mass required in galactic halos or the cosmological background [^. She argued that, 
in the standard non-inflationary scenario, this would happen rather naturally if the 
PBHs formed from scale-invariant fluctuations (n=l) with those of lO^^g having the 
density required to contribute appreciably to 100 MeV cosmic rays. However, the 
relic density would be much reduced if there were an inflationary period. Barrow 
et al. have studied this problem in more detail, calculating the constraints on the 
fraction of the Universe going into PBHs in order that their relics do not have more 
than the critical density |^|. They considered the situation in which PBHs form as 



a result of a flrst-order phase transition induced by bubble collisions at the end of 
an extended inflationary period. We go beyond these calculations in several respects. 
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Firstly, we assume that the PBHs form directly from the inflation-induced density 
fluctuations rather than from bubble collisions. Secondly, we allow for the possibility 
that the equation of state may go soft in some period between the end of inflation 
and reheating, as first suggested by Khlopov et al. . 

The plan of this paper is as follows. In Section 2 we discuss why PBHs may 
form after inflation and why stable relics may be left over from the final stages of 
their evaporation. In Section 3 we derive limits on the fraction of the Universe which 
goes into PBHs of mass M which leave relics, both for the standard case in which 
the Universe is radiation-dominated after inflation and for the case in which there is 
an intermediate dust phase before reheating. In both cases we infer corresponding 
upper limits on the amplitude 5{M). In Section 4 we combine these limits with the 
observed quadrupole anisotropy to infer upper limits on n. If there were no Planck 
mass relics, the strongest upper limit would be associated with the photodissociation 
of primordial deuterium by photons emitted from 10^°g PBHs and this would 
imply n < 1.5 |T^. If there are Planck mass relics and there is no dust phase, the 



limit is strengthened to n < 1.4 for a reheat temperature of order lO^^GeV. If there is 
an extended dust phase after inflation, the limit on n depends crucially on its duration 
and could be as strong as n < 1.3. In Section 5 we draw some general conclusions. 



2 The Formation of PBHs and their Relics 

2.1 The Fraction of the Universe going into PBHs 

We are interested in the situation where PBHs form from the density perturbations 
induced by quantum vacuum fluctuations. We assume spherically symmetric, Gaus- 
sian fluctuations with an rms amplitude 5{M) and a background equation of state 
p = ""fp with < 7 < 1. When an overdense region stops expanding, it must have 
a size of at least times the horizon size in order to collapse against the pressure 
and this requires that 6{M) exceed 7 [^. Thus the probability of a region of mass 
M forming a PBH is 

/3„(M)«i(Af)exp(-^) (2.1) 

and this gives the fraction of the Universe expected to go into PBHs on that scale. 
The mass of a PBH forming at time t must be at least 7^^^^ times the horizon mass, 
so 

M = r^-^Mpi, (2.2) 

tpi 

where rj ~ 7^/^. The value of 7 will usually be 1/3 in the early Universe, corresponding 
to a radiation equation of state. However, one may have 7 = if the Universe ever 
passes through a dust-like phase and, in this case, both Eqs. (|2.1|) and ( p.2|) are 
inapplicable During a dust era, the fraction of the Universe going into PBHs 
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just depends on the probability that regions will be sufficiently spherically symmetric 
to collapse within their Schwarzschild radius and this can be shown to be 



(3{M) ^ 2 X 10-^{My^/^. (2.3) 

If this applies for some period ti < t < t2, then PBHs are expected to form in the 
mass range Mi < M < Mmax, where Mi is given in terms of ti by Eq. ( p.2|) and Mmax 
is the mass of a configuration that just detaches itself from the universal expansion 
at ^2- The latter is calculated as follows. If a region has a fluctuation 6{M) at the 
time tH{M) when it reenters the horizon, then it binds at a time ts = tnS^^^^ with 
a size Rb = ctnS'^. The horizon size at is ctu^'^^'^^ so Rb is 5^^"^ times this and 
the maximal mass is given implicitly by 

M^ax = [S{M^..)f^ Mpi. (2.4) 

The last two terms just give the horizon mass at ^2- In order to determine Mmax 
explicitly, one needs to know the form of 6{M). 

The possibility of a soft equation of state is particularly relevant to the present 
work because, in some circumstances, this may occur during the reheating phase at the 



end of inflation ||30[ and this is precisely the period in which most PBHs are expected 
to form. If inflation ends by means of a second-order phase transition, the scalar field 
oscillates in the potential minimum from the time ti = and then reheating occurs 
as friction generated by the coupling of the scalar field to other matter fields turns the 
kinetic energy of motion into background radiation. The reheating is completed on a 
time scale determined by the decay width F of the scalar field. Providing F"^ <^ ti, 
the reheating time can be taken to be ti, so we can assume that the equation of state is 
hard (7 = 1/3) throughout the period after inflation. In this case, most of the PBHs 
will form at the epoch ti and have the mass Mi = ?7(ti/tpi)Mpi indicated by Eq. 
(p.2[). This is because Eq. ( |2.1| ) implies that, for a blue spectrum, /3o(M) decreases 
exponentially for M > Mi, whereas PBHs with M < Mi (which form before ti) are 
diluted by inflation. Thus we can regard the PBHs as effectively having a 5-function 



mass spectrum, as in the Barrow et al. bubble collision scenario with a mass 
Ml. For F~^ ^ ti, the scalar field undergoes coherent oscillations and the equation 
of state may be soft (7 = 0) from ti until reheating is completed at the time ^2 = E"^, 
when the scalar field decays rapidly into relativistic particles [Q. In this case, the 
PBHs would have an extended mass spectrum, going from Mi to the mass indicated 
by Eq. dU). 



The quantum evaporation of PBHs by thermal emission leads to numerous upper 
limits on the fraction of the Universe going into PBHs at a given time. In the case of 
a radiation equation of state, upper limits on /3o(M) in the range lO^^g < M < 10^ ''g 



have been summarized in Refs. ||T9|, ^ and are shown in Figure (la). The 



constraints on the probability of PBH formation are modified if there is a dust phase 
after inflation (between ti and ^2) because the ratio of PBH density to radiation 
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density no longer increases in this period. The constraints on the fraction I3{M) of 
the Universe going into PBHs during the dust phase are related to [3q{M) via the 



equation 33 



„(M) ^ >UM)rr {^f {^f'\ (2,5) 

This allows the limits on PBH formation during the dust phase to be calculated 
from the constraints which apply if there is no dust phase. The limits on (3{M) are 
indicated in Figures (2a), (3a) and (4a), corresponding to three different choices of 
t2- Comparison of these observational limits with the predicted values of (5o{M) and 
/3(M), given by Eqs. (|2.1| ) and (|2^ ), then leads to constraints on the form of 5{M), 
as shown in Figures (lb), (2b) (3b) and (4b). The constraints on /3o(M), [3{M) and 
5{M) for M > 10^''g are associated with non-evaporating PBHs and correspond to 
the requirement that they have less than the critical density. We also show the limit 
on 5{M) associated with the lack of spectral distortions in the CMB implied by the 
COBE/DMR results |3^; this limit is now stronger than indicated in Ref. |I9| , 



as pointed out independently by Hu et al. [27|. The constraints for M < 10 g are 



associated with Planck mass relics and are the focus of the rest of the paper. 
2.2 Formation of Relics 

The usual assumption is that evaporation proceeds until the PBH vanishes completely 
p8| but there are various arguments against this. For example, the Uncertainty 



Principle implies that the uncertainty in the mass M of a Schwarzschild black hole 
of radius rs = 2GM/c^ must satisfy cAMArs > fi/2 and a black hole evaporating 
below the Planck mass Mpi = {hc/Gy^"^ would violate this ||39|. Also the expected 



energy of the emitted particle exceeds the rest mass of the black hole for M < Mpi. 
For these reasons Zel'dovich proposed that black holes smaller than the Planck mass 
should be associated with stable elementary particles Another argument is 



that quadratic curvature corrections to the gravitational Lagrangian would change 
the Hawking temperature to the form Tbh = kiM~^ — knM'"^, where ki and kn are 
constants, and this becomes zero (implying that evaporation stops) once the mass 
gets down to 

Mrei=f^j Mpi. (2.6) 

Although we do not know what form quantum gravity corrections should take 
as one approaches the Planck mass, it would be surprising if there were none at all. 



The formation of relics is also related to the paradox of information loss 41 
The evaporation of a black hole involves an initially pure quantum state evolving 
into a mixed one and the basic principle of unitarity is thereby violated. To avoid 
such a conclusion, one must either suppose that the information is contained in the 
evaporated particles or that the evaporation terminates when the black hole 
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reaches the Planck mass, the information being stored in a stable or very long-lived 



relic 43 



The above arguments are too vague to be very convincing, but more specific 
arguments for stable relics have been given and we now review some of these. 

Bowick et al. suggest that stable relics may form because black holes can carry 
axionic charge The axionic field does not gravitate and this means that a black 



hole may have arbitrarily large charge without forming a naked singularity. On the 
other hand, it does have a non-zero potential and this would be relevant for string 
interactions (cf. the Aharonov-Bohm effect). This suggests that axion charge will 
become important once a black hole has evaporated down to of order the Planck mass 
since string effects will then apply. Causality and energy conservation presumably 
limit the amount of axion charge that a black hole can radiate within the age of the 
universe, so Bowick et al. infer that it cannot evaporate completely. 

Coleman et al. also argue for a minimum black hole mass on the basis that a black 
hole can have quantum hair which affects its temperature even though it has no clas- 
sical effects They focus on hair associated with a gauge symmetry. Although 



a hole with Zj^ electric charge has no classical hair, quantum effects (associated with 
virtual strings that wrap around the event horizon) generate a non- vanishing electric 
field which decays exponentially with distance from the hole. In fact, this has negli- 
gible effect on the temperature and much more important are the effects of magnetic 
charge. This possibility arises if one considers SU{N)/Zn gauge theories which per- 
mit monopoles since a black hole can then have classical magnetic hair. In this case, 
the Hawking process stops when the temperature reaches zero and this occurs at a 
mass 

Mrei = [n{N - n)/2NY/^e-^Mpi ^ lO^^Mpi, (2.7) 

where n = 1,2, . . . , N — 1 is the Zn charge and e is the coupling constant. These 
relics can still decay to elementary monopoles if the monopoles are light enough. 
Indeed the elementary monopoles may themselves be n = 1 black holes, in which case 
the extreme magnetically charged black holes for other values of n are kinematically 
forbidden to decay to lighter objects of the same charge |Q. Thus the relics are 
stable. 

Gibbons and Maeda consider scale-invariant theories in which Maxwell fields and 
antisymmetric tensor fields are coupled to gravity with a dilaton WTl. The solutions 



in four dimensions are of the Reissner-Nordstrom type, in which the black hole tem- 
perature goes to zero when the mass equals the charge. For a general dimensionality 
D, they find that the thermodynamic behaviour of an electrically and/or magneti- 
cally charged black hole depends on g, the coupling of the dilaton to the Maxwell 
fields. For g > \/D — 3 the black hole evolves to a naked singularity with infinite 
temperature as in the standard picture. For g < V-D — 3 a singular zero-temperature 
endpoint is reached with a mass dependent on the scalar charge and of order g~^Mpi. 
For g = V-D — 3 the hole evolves to a finite temperature final state. It is interesting 
that string theory predicts g = 1 which corresponds to D = 4 in this case. 
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Torii and Maeda consider a theory that couples a scalar dilaton field to a Yang- 



Mills field . It arises as the 4-dimensional effective theory corresponding to higher 
dimensional unified theories. For coloured black holes with zero dilaton coupling they 
find the effective Yang-Mills charge increases as the mass decreases. The temperature 
also increases at low mass and becomes infinite at a mass of O.SSMpi/^fc, where Qc 
is the coupling constant for the Yang-Mills field. At this point the event horizon 
disappears leading to a stable particle-like solution. Similar behaviour occurs for a 
sufficiently small dilaton coupling. They also consider Skyrme black holes and find a 
similar effect. The evaporation shrinks the event horizon until it disappears, leaving a 
non-Abelian particle-like solution. The Skyrmion mass is determined by the coupling 
constants and is of order fsQc, where fg relates to the mass of the Yang-Mills field 



Callan et al. examine field equations which arise from low energy string theory 
50| . The action includes terms quadratic in the curvature tensor, in addition to the 



usual Einstein action, so Eq. (|2.6|) is applicable. Solving the equations of motion 
by perturbation expansion, they find that close to the black hole the dilaton field 
decreases and the string interactions become weaker. The stringy black hole temper- 
ature is unchanged for D = 4, but it is lower than the usual Hawking temperature for 
> 5 by a factor proportional to the inverse string tension. As such black holes lose 
mass, their temperature reaches a maximum before falling to zero at a finite mass. 

Myers and Simon study black hole thermodynamics in second-order theories in 
higher dimensions Here the Lagrangian is the sum of k combinations of the 



Riemann invariants that give second-order field equations. When D = 2k + 1 = 5, 
they find that the temperature can vanish at finite mass. Although the horizon also 
vanishes at this point, it requires an infinite amount of time to evaporate to this limit 
and a naked singularity does not form. Whitt also considers second-order gravity 
theories and obtains similar results 



3 Observational Constraints from PBH Relics 

In this section we will derive the constraints which can be placed on the fraction of 
the Universe going into PBHs of mass M in order to avoid their relics having more 
than the critical density. We will assume that the relics have a mass nMpi, where 
the parameter k is in the range 1 - 10^ for the scenarios discussed above. We note 
at the outset that one can place no constraints on PBHs which form before the end 
of infiation at ti since they will have been diluted away. Thus there are no limits for 
PBHs with mass below M < Mi = ri{ti/tpi)Mpi. Although we do not know the value 
of ti a priori, we can place an upper limit T^ax on the reheat temperature and thus 
Ti from the observed CMB quadrupole anisotropy. Thus there are no constraints in 
the mass range 

M < M^^ = 7]Mpy{T^^jTpy)-^ (3.1) 
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If one assumes that the quadrupole anisotropy is due entirely to tensor (gravitational 
wave) fluctuations, and that the spectral index n does not deviate significantly from 
unity, the expansion rate of the Universe during inflation cannot exceed H = 2.9 x 
10~^Mpi p3| , This leads to a maximum reheat temperature of 1.6 x 10^^ GeV 
and a minimum mass Mmin ~ 2g. A more careful determination of My^m, allowing 
for the contribution of the scalar fluctuations to the quadrupole anisotropy, is given 
in Section 4 [c.f. Eq. (^)]. 

3.1 Constraints for a Hard Equation of State 

Let us first consider the F^^ <^ ti situation, in which the early Universe is radiation- 
dominated immediately after inflation. The form of the relic constraint depends 
crucially on whether the PBHs dominate the density before they evaporate at |^ 

tevap(M) = (M/Mpi)3tpi. (3.2) 

Since the ratio of the PBH density to radiation density increases as t^^^, the condition 
for the radiation to dominate at evaporation is 



where t is the time at which PBHs of mass M form and is given by Eq. ( p. 21) . This 
assumes that there are no PBHs in a different mass range; we are essentially taking 
the PBHs to have a 5-function mass spectrum. In this case, the ratio of the relic 
density to the critical density at the present epoch is 

"..^f.^)f^)f^)". (3.4) 



I- \ M J \t 

where tgq = ^o^rad is ff^^ time at which the matter and radiation densities are equal 
and to = 6.5/i~^Gyr is the age of the Universe for VLq = 1. Here Q^ad = 2 x 10~^/i~^ 
is the current radiation density parameter and h is the Hubble parameter in units of 
100 km s~^ Mpc^^. Using Eq. (p^.2|), the constraint firei < 1 becomes 



^<,r"V'(^)'%^'^(^)''^:o-v-.-v(-^)'^ ,3.) 



The PBHs dominate the density at tevap if condition (|3.3| ) is not satisfied. In this 
case, most of the background photons derive from the PBHs, so Eq. ( p.4| ) is replaced 
by 

Mpi \ (t^\ ^ fA£\ ft^\ V2 
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and the constraint firei < 1 becomes 



M > K^/^ ( ^) f]J/]°Mpi = 10'\'/'h-^/'Mp,. (3.7) 
\tp\J 

Note that this hmiting mass also corresponds to the intersect of conditions ( p.3|) and 
( p.5|) . Eqs. (|3.1j ), (|3.5| ) and (p.7|) define the shaded fines on tfie left-fiand-side of 



Figure (la). Tfiese equations were impficitly derived by Barrow et ah p9|, afifiougfi 
tfiey expressed tfiem in terms of T ratfier tfian M. Tfiere is also a mistake in tfieir 
Eq. (5.20) and tfiis leads to an error in tfieir Eq. (5.31) wfiicfi is equivalent to our 
Eq. (P). 

Tfie relics fiave tfie critical density only on tfie boundary specified by Eqs. ( |3.5| ) 
and ( |3.7| ). It is interesting tfiat any value of P above 10~^^ will suffice to provide 
tfie critical density if M is fine-tuned to fiave tfie value given by Eq. ( p.7|) . Recall 
tfiat most of tfie PBHs actually form at tfie end of inflation, i.e. witfi tfie mass 



given by Eq. (3J.), so tfiis corresponds to fine-tuning tfie refieat temperature to 
T = 10~^r7^/^K~^/^Tpi ^ lO^^GeV. Altfiougfi one migfit regard tfiis as unlikely, all tfie 
present-day radiation originates from PBH evaporations in tfiis situation, so in tfiis 
case one migfit regard firad as a free parameter determined by tfie refieat temperature. 

Two otfier constraints are sfiown on tfie left of Figure (la). Tfie stronger limit 
(sfiown dotted) comes from tfie fact tfiat tfie observed baryon asymmetry of tfie 
Universe could not be generated below tfie critical temperature Tmin = O(10^)GeV 
associated witfi tfie electroweak pfiase transition. PBHs tfiat form witfi a mass in 
excess of 

Mpi(T^i,/Tpi)-2/3 = lO^g (3.8) 

evaporate after tfiis, so, if sucfi fioles come to dominate tfie Universe before evapora- 
tion, tfie evaporated radiation is not sufficiently fiot to allow baryogenesis to proceed. 
Tfiis implies tfiat condition ( p.3|) must apply for M > lO^g and tfiis corresponds to 
tfie dotted line in Figure (la). Coincidentally, tfie masses given by Eqs. ( |3.7| ) and 
( p.8| ) are nearly tfie same, so tfiis would exclude tfie critical density condition associ- 
ated witfi Eq. (|3.7|) . However, tfiis limit is not completely secure since tfiere may be 
otfier mecfianisms (including black fiole evaporations tfiemselves [^) for generating 
baryon-asymmetry. A weaker but more reliable constraint corresponds to tfie require- 
ment tfiat evaporating PBHs do not generate a pfioton-to-baryon ratio exceeding tfie 
current value 5*0 = 10^ One can sfiow tfiat tfiis just corresponds to tfie condition 



^0 <109f-^y' (M<10"g), (3.9) 



'PI 



wfiere tfie upper limit on M arises because tfie PBHs must evaporate early enougfi 
for tfie pfiotons to be tfiermalized. Tfiis is weaker tfian condition ( |3.3| ) by tfie factor 
5*0 and is labelled 'entropy' in Figure (la). 
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In order to convert these constraints on (3o{M) into constraints on the horizon- 
scale density fluctuations S{M), we use Eq. (|2.1| ) to obtain 



S{M) < 0.16 [logio5(M) - logio/5o(M)] 



-1/2 



(3.10) 



where the value of (3o at a given M is determined by the boundary of the shaded 
regions in Figure (la). The relic constraint then becomes 



6{M) < 0.16 



0.13 



27-logio(/«r7i/') + 1.51ogio 



/ M 



+ logio5(M) 



-1/2 



17 - l06 



•10 



Ml 



Pi 



for M < 10"/t2/5/i-4/5Mpi, 



(3.11) 



where we assume k, = t] = 1 and take \ogiQ6{M) ^ —1.5 in the second expression. 
The relic limit on 6{M) is shown on the left of Figure (lb). 



Figures la &: lb 



3.2 Constraints for a Soft Equation of State 

We now consider the case in which the Universe has a dust era between ti and t2- 
We first assume that this ends before the PBHs evaporate, so that 



(3.12) 



Mpi \tpiJ 

from Eq. ( p.2|) . The condition that the radiation dominates the density when the 
PBHs evaporate now becomes 



1 — /5 \tcvap J Vtpi/ VMpi. 

with no dependence on the time at which the PBHs form and the condition i^rei < 1 
corresponds to the upper limit 



^ < 10-'^ K-'Ji' . (3.14) 

1-/3 \MpJ \tpJ 

Eqs. ( p.l3| ) and ( p. 141) just come from Eqs. (|3l3| ) and (|3.5| ) together with Eq. (|2.5|) . 
If the PBHs start to dominate the density before tevap, then Eqs. ( p.6|) and ( ^.71) still 
pertain. 

Next we assume that the dust era ends after tevap, so that Eq. ( p.l2|) is violated 
and the density of the post-inflationary Universe is never dominated by radiation 
before ^2- The current relic-to-radiation ratio therefore becomes 

't \2/31-^ 

''evap \ 



M J \ t 



(3.15) 
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where the two terms in square brackets give the contributions to the radiation density 
from the reheating and the PBH evaporations. Since the latter can be neglected, the 
condition firei < 1 is again given by Eq. ( |3.14|) . 



In order to determine the general constraints on (3{M) and 6{M) one needs to 
specify the epoch t2 at which the dust phase ends. Eq. ( p.4| ) then determines the 
mass range of PBHs forming during the dust era. For Mmin ^ M < M^ax the 
constraint on /3(M) is given by Eq. (|2.5D, where the value of I3q{M) is determined by 
the shaded boundary in Figure (la). Eq. (p.3|) then implies the upper limit 



, / X 1/13 / ]\,f N -1/13 

5{M) < l.mmf"' (^) V'"' (3.16) 

and, in particular, the relic limit given by Eq. ( |3.5| ) implies 

. f.,T s 2/13 / / X 1/13 

6{M) < 1.3 X IQ-^K-^/'^h'/^^ (^] ( ^ ) . (3.17) 

\M-p\J \tpi/ 

The relic limits on (3{M) and 5(M), together with the previously known constraints, 
are shown in Figures (2), (3) and (4). Figure (2) applies if the dust era ends with 
the formation of lO^g PBHs (^2 = 10~^^s), corresponding to the mass given by Eq. 
(p.7|). Figure (3) applies if it extends until 10~^^s, which is just after the 10^°g PBHs 
form. Figure (4) applies if it extends until 10~^^s, which is just after the lO^^g PBHs 
form. [Eq. ( P^ implies that the condition for this is t2 > 10"^^5(10^^g)~^/^s; this 
exceeds the usual time of 10~^'^s because the PBHs are now smaller than the particle 
horizon at formation.] Eqs. ( |3.1D , ( |3.7D and ( |3.14| ) define the shaded lines on the 
left-hand-side of Figures (2a), (3a) and (4a). Qualitatively, the effect of softening is 
to bring down the constraints relative to the hard case. For low values of t2, only 
the relic constraint is brought down, as indicated in Figures (2). As ^2 increases, the 
deuterium and gamma-ray limits are also brought down, as indicated in Figures (3) 
and (4). 

Figures 2, 3 & 4 



4 Results for a Constant Spectral Index 
4.1 Normalization to the Quadrupole 

So far the analysis has been completely general and has made no assumptions about 
the form of the power spectrum of the density fluctuations. In this section we shall 
consider the special case in which the spectral index n is constant and consider the 
corresponding upper limit on its value. For a general infiaton potential, n will not be 
exactly constant but we show in the Appendix that this is always a good approxima- 
tion in the hybrid inflationary scenario. The strongest upper limit on n is derived by 
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normalizing at the COBE/DMR quadrupole scale and finding the steepest straight 
line which avoids all the shaded areas in Figures (lb), (2b), (3b) and (4b). In case (1), 
the strongest constraint is associated with the relics from the PBHs of mass Mmin. 
This means that relics can have the critical density but that PBHs cannot contribute 
appreciably to lOOMeV cosmic rays. In case (2), the constraint on n associated with 
the relics is even stronger because only the rehc dip is brought down by the dust era. 
However, in case (3), the deuterium limit gives the strongest constraint on n because 
the 10^°g dip has also been pulled down. In case (4) the gamma-ray limit gives the 
strongest constraint because the lO^^g dip has been pulled down, so PBHs can gen- 
erate the cosmic rays but not a critical density of relics. Note that it is not possible 
to produce both the cosmic rays and a critical relic density, i.e. there is no situation 
in which a single straight line passes through both the Mmin and lO^^g points. Either 
the lO^^g PBHs form during the dust era, in which case the gamma-ray limit implies 
the relics are unimportant, or they form after the dust era, in which case the relic 
limit implies the gamma-ray limits are unimportant. 

In order to determine the constraint on n, we need to normalize the density spec- 
trum at the quadrupole scale. To do this, we must first convert the observed tem- 
perature fluctuation into a corresponding density fluctuation. Inflation also produces 
primordial tensor fluctuations (gravity waves) and these can be significant in some 
cases. If no reionization occurs, the surface of last scattering is located at a redshift 
zi^s ~ 1100 for i7o = 1 and the angle subtended by the horizon at that redshift is ap- 

— 1/2 

proximately 6 ^ z^^^ ~ 2°. Once gravity waves reenter the horizon, they redshift as 
relativistic matter and soon become negligible. However, experiments probing angu- 
lar scales larger than 2° measure superhorizon-sized perturbations and are therefore 
sensitive to gravitational wave effects. 

The CMB anisotropics can be expanded into spherical harmonics with coefficients 
aimi'x) which are stochastic random variables. The observed multipoles detected from 
a single point in space are defined in terms of these coefficients: Qf = ^ Y.m=-i k«mP 
Il57| . A given inflationary model predicts values for the averaged quantities 6f = 
^(2Z -|- l)(|aim(x)p), where the average is taken over all observer points. If the 
aimi'x.) are uncorrelated, the expected contributions from the scalar (6f [S*]) and tensor 
(6f [T]) fluctuations add in quadrature, so 6;^ = QflS] + O^^]- 



We follow the notation of Ref. [0 and define the scalar and tensor contributions 
to 62 as 

s.e5|s, = 5%M), T.eri = 5%EQ. (4.1) 

4:71 Atc 
respectively. The observable quantities S and T have been calculated numerically 
in terms of the corresponding values of the potential and its derivatives |T^. The 
relevant expressions for the present work are 

-1.65 fl + 0.20^')T (4.2) 
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Ml, 



3.14 



2^3/2 



(4.3) 



where a subscript C refers to the quadrupole scale. Substituting these expressions 
into Eq. (|1.1| ) yields the magnitude of the density fluctuation at the quadruole: 



5c ~ 0.675^/2 (1 + 0.3 



(4.4) 



4.2 Constraints on the Spectral Index for no Dust Phase 

If the Universe has the critical density, the comoving rest mass within a sphere of 
radius A/2 is M = 1.5 x lO^^M0/i^(A/Mpc)^. The quadrupole corresponds to a scale 
Ac ~ 6000/i~^ Mpc and has an associated mass Mc = 10^^h~^g. For a power law 
spectrum we require 5c{M / Mc)^^'"^^^^ be less than the value of 5{M) indicated in 
Figures (lb)-(4b). Taking logarithms and substituting Eq. ( [4. 41 ) then yields an upper 
limit on n in terms of the 5* and T: 



n- 1< 6 



log 



10 



M 



-1 r 



0.17 + log 



10 



- log 



10 



T 

1 + 0.3^ 



. (4,5) 



The last logarithmic term is a correction accounting for the gravitational wave contri- 
bution to the observed quadrupole anisotropy. Since \ogiQ{5 / S^^"^) > 4, we conclude 
that this correction term is negligible. Furthermore, although the COBE/DMR ex- 
periment measures the quantity S + T, it is easily seen that limit (|4.5| ) is not altered 
significantly if it is assumed that the entire COBE/DMR signal is due to the scalars. 
This shows that a separate determination of the gravitational wave contribution is 
not required in order to derive a limit on n. 

A second feature of this limit is its relative insensitivity to cosmic variance. Cosmic 
variance arises because the fluctuations predicted by theory are stochastic in nature 
and have a Gaussian probability distribution. A set of observations can only measure 
a finite number of realizations of this distribution and this is never sufficient to spec- 
ify it completely. Thus there always exists an intrinsic uncertainty in the analysis. 
Although the effect is most significant at the quadrupole scale and can be problematic 
when normalizing large-scale structure observations, it should not significantly affect 
the PBH limit since this spans such a large range of scales. 

It only remains to substitute the numerical value of S^^'^ into Eq. ( [4.5| ). A max- 
imum likelihood estimation made by Seljak and Bertschinger ||5^ with the first-year 
COBE/DMR data implies a central value of 



5^/2 = 5.8 X io-6eO-46(i-") 



(4.6) 



if the blackbody temperature of the CMB radiation is T^ad = 2.736K. Smoot et al. 



and Bennett et al. [12| arrive at similar expressions, although we emphasize that 
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the precise form of the relationship between S and n is not important for the PBH 
constraints. Hence, substituting Eq. (|4.6|) into Eq. ([4.5|) imphes 



An expression for M^[^ may be derived in terms of quantities that are in principle 
observable. During reheating the potential energy of the scalar field driving inflation is 
converted into relativistic particles with an energy density prad = where e is 

the number of relativistic degrees of freedom. Typically e = (9(10^) and e ~ 160 in the 
simplest SU(5) GUT model. An observable upper limit on the reheat temperature 
may be calculated by assuming that the available energy is given by Vc- For an 
efficient reheating process this gives 



30.. 



Ti < ( —Vc 



, ^ + 0.2-)T 

e V S 



Tpi, (4.8) 



where the second expression uses the lowest order terms in Eq. ( |4.2| ). The smallest 
black holes formed after inflation must therefore have a mass exceeding 

M^i, = 1.7 X 10-6 i^^y (l - 0.l|) g. (4.9) 

In principle one requires a knowledge of the tensor contribution to determine Mmin 
but such information is not currently available. 

In the efficient reheating case, the strongest upper bound on the spectral index 
n depends on the value of Mmin or equivalently ti. For Mmin < 10^ g, the strongest 
limit is associated with the relics and the relevant mass is Mmin itself. Substituting 
Eq. (|3?TT1) into Eq. (^) then gives 

n-l< , [4.5 - ^log.o fl2 - log.o (MlnAVi . (4.10) 

56 - logio(Mmm/g) L 2 V \ g //. 

Thus the limit on the spectral index is n < 1.4 for Mmin ~ Ig? corresponding to 
Ti = O(10i6)GeV, and n < 1.5 for M^in ~ lO'^g, corresponding to Ti = O{10^^)GeV. 
The limits on n for other ranges of Mmin are shown by the upper shaded curve in 
Figure 5. Here we are regarding n as a continuous function of Ti (or equivalently ti). 
The closure density constraint applies for Mmin > 10^'^g (ti/tpi > 10^^), the gamma- 
ray constraint for lO^^g < Mmin < lO^^g (10^*^ < ti/tpi < 10^^) and the deuterium 
constraint for 10^°g < Mmin < lO^^g (10^^ < ti/tpi < 10^^). Since the helium 
and entropy constraints are associated with higher values of (5h than the deuterium 
constraint, the limit on n does not improve until Mmin gets down to lO^g when the 
relic constraint ( [4.10| ) becomes relevant. This is why the part of the upper curve in 
Figure (5) between "relics" and "deuterium" is flat, with a discontinuity at lO^g. Note 
that the CMB distortion gives a better limit on n than the closure density limit for 
Mmin > lO^^S' (or ti > lO^^^s) and this is why the limit flattens off at n = 1.76. Hu 
et al [0 have calculated this limit somewhat more carefully and flnd that it flattens 
off at n=1.54, as indicated by the dotted line in Figure 5. 
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4.3 Constraints for an Early Dust Phase 

If there is a dust phase after inflation, the constraints on the index n depend on both 
the value of ti and ^2, with the latter itself being determined by the decay width V. 



We first assume that ti is fixed at the value associated with Eq. ( [4.8| ) and allow F to 
vary. An upper limit on F follows from the fact that the reheat temperature at the 
end of the dust phase is Trh ~ (rtpi)^''^Mpi. Therefore if efficient conversion of the 
vacuum energy to relativistic particles occurs, the assumption that the COBE/T^MK 
quadrupole signal is due entirely to gravitational waves implies Ftpi < 10~^ as shown 
by the shaded line on the left of Figure 5. The more general limit is Ftpi < (Ti/Tpi)^ 
where Ti is given by Eq. ( [4.81 ). A lower limit on Trh follows from the requirement that 
baryogenesis must proceed after reheating. If the lowest temperature for which the 
observed baryon asymmetry may be generated is the electroweak scale, (9(10^) GeV, 
we require Ftpi > 10"^*^. As discussed in Section 3.1, there is some uncertainty in the 
argument, so we show this limit by a broken limit in Figure 5. For intermediate values 
of F, one can merely exclude certain areas in the (F, n) plane as we now demonstrate. 

If we normalize on the COBE/J^MR quadrupole scale, Mq-, the rms amplitude 
on a smaller scale M is 5{M) = 5c{M /McY^'''^/^ where 5c ~ 3.8 x 10"^ from Eqs. 
( [4.4|) and ( [4.6|) . Substituting this expression into Eq. ( p. 41) implies that the maximum 
mass of a PBH formed during the dust phase is 

. /2\4/(n+3) („_,)/(„+3) 

It follows that PBHs with mass M are formed during the dust phase only if 

logioHpi < -23.6 + 15.5n - log^o (^) • (4.12) 



Another upper limit on Ftpi follows from Eq. ( p.l6|) with ^2 identified with F ^: 



logioHpi < 208 - 134n + ^(13r2 - 19)logio (^) + 21ogio/3o(M) + log^or/, (4.13) 

where the value of [3q{M) is indicated in Figure (la). For a given value of F, Eqs. 
( [4.12| ) and ([4.13|) can also be interpreted as giving lower and upper limits on n. 



The strongest limit on n or F is associated with the relic constraint for M^ax < 
10^°g, the deuterium constraint for 10^°g < Mmax < lO^^g and the gamma-ray con- 
straint for Mmax > lO^^g. Eq. ( [4.11|) implies that the relic constraint applies for 

logioFtpi > -35 + 12n (4.14) 

and Eqs. (U) and (|1|) with M = M^in ~ Ig then imply 

logioFtpi < 153 - 125n. (4.15) 
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Similarly the deuterium constraint, /3o(10^°g) < 10 applies for 

- 35 + I2n > logio(npi) > -38.6 + 10.5n (4.16) 

and Eq. (^A3|) with M = 10^°g gives 

logio(rtpi) < 118 - lOln. (4.17) 

The gamma-ray constraint, /?o(10^^g) < 10^^^, applies for 

logio(rtpi) < -38.6 + 10.5n (4.18) 

and Eq. (^J3D with M = lO^^g gives 

logio(rtpi) < 93-91n. (4.19) 

These constraints on (r,n) are shown by the lower shaded line in Figure (5). The 
relic constraint applies for 10~® > Ftpi > 10^^'', the deuterium constraint for 10^^'' > 
rtpi > 10^^^ and the gamma-ray limit for 10~^^ > Ftpi > 10~^°. We therefore arrive 
at an upper limit of n = 1.4 if PBH form relics and there is an extended dust phase 
immediately before the standard radiation-epoch. This limit is independent of the 
form of the dark matter, the bias parameter, cosmic variance and any gravitational 
wave contribution to the COBE/DMR signal. It is clear from Figure (5) that the range 
1.3 < n < 1.4 is astrophysically very interesting. For Trh < lO^GeV the origin of the 
observed gamma-ray and cosmic ray spectra can in principle be explained in terms 
of evaporating lO^^g PBHs formed during the early dust phase. This corresponds 
to the limit (^4.19 ) becoming an equality. For Trh > 10^'^GeV, however, the Planck 



mass relics of Ig PBHs may be a natural candidate for the cold dark matter in the 
universe. Note that if ti is allowed to increase associated with Eq ( [4.8| ), then the 
limits on (n, F) are somewhat weaker than indicated by the lower curve. The lower 
curve therefore gives the most stringent limit. 

Figure 5 



5 Conclusions and Discussion 

In this paper we have derived upper limits on the fraction of the Universe going into 
PBHs at various mass scales and thereby inferred upper limits on the rms amplitude of 
the initial density perturbations. In the inflationary scenario the spectrum of density 
fluctuations arising from quantum fluctuations in the vacuum is uniquely determined 
by the functional form of the self-interaction potential of the inflaton fleld, so it is 
possible to constrain the inflaton potential with the PBH constraints. 

One of our main purposes has been to explore the consequences of the suggestion 
that evaporating PBHs may leave behind stable relics with masses of order the Planck 
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mass. The requirement that such rehcs have less than the critical energy density leads 
to constraints which are stronger than those previously derived for masses smaller 
than lO^g. We have also investigated to what extent the PBH constraints are altered if 
there was an early dust phase in the history of the Universe. Generically inflation ends 
via a phase transition which can be either first- or second-order. In the latter case, the 
scalar field rolls down to the global minimum of its potential and begins to undergo 
coherent oscillations. If the decay of the field to relativistic particles is rapid relative to 
the expansion rate of the Universe, reheating occurs almost instantaneously. However 
these oscillations can be prolonged if the coupling constants are sufficiently small 
and, in this case, the Universe becomes dust like for some period before the standard 
radiation epoch is recovered. This is very important for the PBH constraints because 
this is precisely the period in which the PBHs are expected to form. The constraints 
now depend crucially on the duration of this dust phase but the upper limits on the 
rms amplitude are always strengthened. 

In order to constrain the spectral index n of the density fluctuations, we have 
normalized the rms amplitude on the COBE/J^MK quadrupole scale and applied the 
above results to the special case in which the spectral index is constant over all the 
scales of interest. This example has recently been investigated within the context 
of the hybrid inflationary scenario. The limit on n depends crucially on the reheat 
temperature Trh- If there is no early dust phase, the relic constraint provides the 
strongest upper limit for Trh > lO^^GeV and this corresponds to PBHs with masses 
smaller than lO^g. The limit is n ^ 1.42 for Trh ~ lO^^GeV and n ^ 1.49 for 
^RH ~ lO^^GeV. The constraint is tightened if there was an early dust phase and 
it was found that the range 1.3 < n < 1.4 has a number of potentially interesting 
astrophysical consequences. In particular relics can provide the critical density for 
Trh > 10^'^GeV and lO^^g PBHs can contribute appreciably to the observed gamma- 
ray and cosmic-ray spectra at lOOMeV for Trh > lO^GeV. It is interesting that the 
current best-fits to the second-year COBE/DMR data appear to favour this range for 
the spectral index, although it must be emphasized that these fits are also consistent 
with a flat (n = 1) spectrum. It is also intriguing that the BATSE detector on the 
Compton Gamma Ray Observatory (CGRO) has observed gamma-ray bursts that 
appear to exhibit the characteristics associated with PBH explosions at the present 
epoch |61|. It remains to be seen, however, whether this range of values for the 
spectral index can be reconciled with current data from large-scale structure. 

A number of assumptions have been made in this work. For simplicity we have 
assumed that the transition from the coherent oscillation phase to the radiation phase 



occurs instantaneously, although this is not necessarily the case [34|. Moreover, during 



the oscillation phase, the scalar field only behaves as a dust fluid on average over one 
cycle. The formation of PBHs in such a case may be different to that assumed in Ref. 
||33| and presumably this would affect the observational limits. Moreover, in deriving 
the constraints on the spectral index, we assumed that n is constant over all relevant 
scales. 
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In conclusion, therefore, the above constraints are useful for a number of reasons. 
Firstly they extend to scales many orders of magnitude smaller than those accessible 
to large-scale structure observations. This imphes that errors in the COBE/DMR de- 
tection due, for example, to cosmic variance are not important when one normalizes 
the spectrum on the scales associated with large angle CMB experiments. Likewise, 
any contribution to the large-scale CMB anisotropy from a primordial spectrum of 
gravitational waves is generally also negligible. Moreover, both the COBE/DMR 
point and the PBH limits are independent of the form of the dark matter in the 
Universe and the constraints on the spectral index are therefore independent of any 
transfer functions and biassing parameters. Consequently, although they are some- 
what weaker than those derived from current large-scale structure observations alone, 
they have the advantage that they are relatively insensitive to any specific choice of 
cosmological model. 
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Appendix 



In the case of hybrid inflation, the constancy of the spectral index is a natural result 
P^ , In this scenario one of the scalar fields is initially at zero and the potential 
takes the asymptotic form 

y(O,0) = \/o + ^mV, (A.l) 

where Vq = M^/4A^ and M is a free parameter. Such models assume Am^ <^ M^g"^ 
and SXm'^Mpi <^ 27rM^, where A and g are coupling constants typically of order unity 



2^ . The inflationary phase ends when reaches 0e = M/g and triggers the rapid 



rolling of the second field. The spectral index and number of e-folds are given by ||T6| 

^ = i + MZ!_r^(_) (A.2) 

(A.3) 




where the quantities on the right hand side of Eq. ( [A.2|) are evaluated when the 
scale first crosses the Hubble radius during inflation. If one requires n > 1, the false 
vacuum term Vq in the potential must dominate the m^0^ contribution |jT9[. This is 
the case near the end of inflation, so the changes in n and N are given by 

An ~ _Pi_^A0 A.4 

Air K) 



and 



AiV ^ 77^zlzA0, (A.5) 



Sir Vo 

respectively. Eliminating A0 and substituting for Vq with = 0e gives 



A ^rS^rA^ (mMpy /saV Am^ 
An^-AA^— (^^-j ^ (A.6) 

and for AN ~ 60 this implies that An ^ 1 over the range of interest. We conclude, 
therefore, that n is effectively constant over the last 60 e-foldings of inflation. This 
may be understood physically by considering the form of the effective potential ( |A.1| ) 
as rolls towards the minimum. As decreases the form of Eq. ( |A.1|) appro aches 
the secant potential that leads to an exactly constant spectral index. 
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Figure 1: (a) The constraints on the fraction (3q{M) of the early Universe going 
into PBHs with mass M if the equation of state is radiation-hke (7 = 1/3). The origin 
of the constraints above lO^'^g is summarized in Ref. |19]. The "entropy" constraint 
arises from the requirement that the PBH evaporations do not generate more than 
the observed photon-to-baryon ratio. There is potentially a stronger but less secure 
constraint in this mass range derived from the assumption that the observed baryon 
asymmetry must be generated above the electroweak scale and this is shown dotted. 
Below lO^g the strongest limit is due to the relics left over from PBH evaporations, 
(b) The corresponding constraints on the rms amplitude of the density fluctuations. If 
PBHs do not leave behind stable relics after evaporation, the strongest upper bound 
on the spectral index is given by the dashed line which joins the COBE/'DMR point 
and the deuterium constraint at lO^^g. This limit applies for reheat temperatures 
above (9(10^)GeV. If relics are formed, the limit is strengthened at higher reheat 
temperatures as indicated. 

Figure 2: (a) The constraints on /3(M), the fraction of the Universe going into 
PBHs of mass M during a post-inflationary dust phase. The dust phase is assumed 
to last until t2 ~ 10~'^^s, long enough for lO^g configurations to detach themselves 
from the cosmological expansion before the radiation era begins, (b) The correspond- 
ing constraints on 5(M), with the lines having the same significance as in Figure 
(lb). Although the limit on the probability of PBH formation is reduced during 
the dust phase, the upper limit on the rms amplitude of the fluctuation at a given 
scale is increased because the fraction of the Universe going into PBHs is no longer 
exponentially damped. 

Figures (3a) and (3b): The same as Figures (2a) and (2b) with the duration of 
the dust phase extended until ^2 ~ 10~^^s to allow 10^°g PBHs to form. 

Figures (4a) and (4b)'. (a) The same as for Figures (2a) and (2b) with the duration 
of the dust phase extended until ^2 ~ lO^^^s to allow lO^^g PBHs to form. 

Figure 5: Illustrating the constraints on the spectral index arising from the over- 
production of primordial black holes, the shaded area being excluded. The lower line 
applies if there is a dust phase immediately after inflation, in which case the ordinate 
is logioT^pi, the upper line if there is no dust phase in which case it is log]^Q(tpi/ti). 
The constraints depend on the reheat temperature Trh ~ 10^^(rtp)^/^GeV, where 
r is the decay width of the scalar field that decays into relativistic particles. The 
n-independent upper and lower limits on the decay width arise from assuming the 
COBE/'DMK detection is due entirely to gravitational waves (shaded line) and from 
requiring that baryogenesis can only proceed above the electroweak scale (dashed 
line). The dotted horizontal line indicates the CMB distortion limit of Hu et al | |27|| . 
For reheat temperatures above Trh ~ 10^'^GeV the most important constraint arises 
from the requirement that any Planck mass relics left over from the final stages of PBH 
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evaporation should have less than the critical density at the present epoch. For lower 
reheat temperatures, more massive PBHs may form and the strongest constraints 
then arise from the photodissociation of deuterium by evaporating 10^°g PBHs, from 
the observed gamma-ray background in the energy range 0.1 — 1 GeV, or from the 
distortions of the CMB. It is clear from this figure that the region 1.3 < n < 1.4 
has a number of interesting astrophysical consequences if there was an extended dust 
phase immediately after the inflationary expansion. 
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